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Which is the flavor of cosmic neutrinos seen by IceCube? 
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We analyze the high-energy neutrino events observed by IceCube, aiming to probe the initial 
flavor of cosmic neutrinos. We study the track-to-shower ratio of the subset with energy above 60 
TeV, where the signal is expected to dominate and show that different production mechanisms give 
rise to different predictions even accounting for the uncertainties due to neutrino oscillations. We 
include for the first time the passing muons observed by IceCube in the analysis. They corroborate 
the hypotheses that cosmic neutrinos have been seen and their flavor matches expectations. 
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I. INTRODUCTION 

The search for High Energy Starting Events (HESE) in 
IceCube detector provided the hrst evidence for a high- 
energy neutrino flux of extraterrestrial origin [THSl- In 
three year of data taking [T], 37 events with deposited 
energies above 30 TeV were observed, relative to an ex¬ 
pected background of 8.4 ± 4.2 cosmic ray muon events 
and 6.6±5.9 atmospheric neutrinos. 

The scientific debate about the origin of these events 
is extremely lively. There is little doubt that cosmic neu¬ 
trinos have been seen, but their origin and propagation 
is not understood. In order to proceed, the flavor com¬ 
position has to be probed. The flavor discrimination is, 
in principle, possible by looking at the topology of the 
events. Most HESE have ‘shower’ topology, that includes 
neutral current (NC) interactions of all neutrino flavors 
and charged current (CC) interactions of Vg and being 
the decay length of the r lepton too short to be resolved 
below ^ IPeV. On the other hand, events with ‘track’ 
topology are produced by CC interactions of v^. Thus, 
the crucial observable quantity is the ratio between track 
and shower events at high energy and it can be used to 
confirm the cosmic origin and/or to discriminate among 
different production scenarios. With this purpose the 
authors of mu recently discussed the observed track-to- 
shower ratio of the IceCube data with energy above 30 
TeV claiming that these data are marginally compatible 
with the hypothesis that neutrinos are of cosmic origin. 
These studies have been influential, setting the case for 
a muon deficit problem in IceCube, see e.g. mm- 

In view of the importance of this issue, we perform 
an independent analysis adding our contribution to the 
discussion. We focus on the subset of events with de¬ 
posited energy above 60 TeV, where the signal is expected 
to dominate. We show that different production mecha¬ 
nisms give rise to distinctive expectations of the track-to- 
shower ratio, even when the uncertainties due to neutrino 
oscillations are included. Also the muon neutrinos pass¬ 
ing through the Earth confirm the existence of an astro- 
physical component and we include for the hrst time this 
information on the analysis. We hnd that present data 


set is well compatible with the hypothesis that cosmic 
neutrinos have been seen, even if the limited statistics 
does not allow yet to discriminate the initial havor. 

II. EXPECTATIONS 

From neutrinos to HESE events. Let us consider 
HESE events with energies between 60 TeV to 3 PeV 
and starting inside IceCube, that are likely to be dom¬ 
inated by the signal due to cosmic neutrinos. The ex¬ 
pected number of events produced by an isotropic flux 
of neutrinos and antineutrinos with havor £ is, 

N = 4Tr T J dE A?e{E)Ae{E) (1) 

where £ = e, r and T is the observation time. The 
effective areas Ai{E) are provided by the IceCube col¬ 
laboration [3] and include the ehects of neutrino cross 
sections, partial neutrino absorption in the earth, detec¬ 
tor efficiency and specihc cuts of the HESE analysis. 

In order to calculate the track-to-shower ratio, we sep¬ 
arate the different contributions to the effective areas, 

A^{E) = AliE) + AliE) (2) 

where Aj^{E) = pt{E) A^[E) is the effective area for 
CC interactions that produce tracks in the detector, while 
A^E) ^ (1 — Pt(A)) A^{E) is the effective area for neu¬ 
tral current (NC) interactions that are instead observed 
as showers. The parameter pt{E) is the probability that 
an observed event (i.e. passing all the cuts in the HESE 
analysis) produced by a muon neutrino with energy E is 
a track event. This quantity is given by, 

’ anc(E)M«C(E) + „cc{E)MC°(E) 

where crcc{E) and crNc(A') are the cross section for CC 
and NC interactions of neutrinos 0 while M^'^[E) and 
are the effective detector mass for CC and NC 
interactions oi [3]. The probability px is mildly de¬ 
pendent on energy and approximately equals 0.8. 
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The number of showers Ng and tracks in the Ice- 
Cube detector can be then calculated according to: 


where U is the neutrino mixing matrix and are the 
flavor fractions at the source given by: 


pE 

Ns=iTTT dE {<^e{E)A,{E)E<^r{E)Ar{E) + 

JO 

+ %{E)[1-pt] A^{E)} 

pE 

Nt=4ttT dE<i>^{E)pTA^{E) (3) 

Jo 

In the above relation, we neglected the small fraction 
of Ur CC-events followed by taus decaying into muons 
which can be potentially observed as tracks. Moreover, 
we introduced an upper integration limit at i? = 3 PeV 
since the HESE analysis only includes events with de¬ 
posited energy below 3 PeV. In principle, the effects of 
the threshold at E^ep = 3 PeV should be implemented 
as a correction of the effective areas. Here, we assume 
that this can be mimicked by a sharp cut in the AflE) 
at the neutrino energy E = 3 PeV. We tested the validity 
of this approach by comparing our predictions with the 
expected numbers of events in Supp. Tab. IV of [T]. We 
obtain good agreement both for the absolute and relative 
numbers of shower and track events. 

Description of cosmic neutrinos. Cosmic neutrinos 
are surely due to non-thermal processes. Thus we ex¬ 
pect that their fluxes averaged over the directions, are 
approximated by a power law distributions up to a max¬ 
imum value that we assume being larger than 3 PeV, 


<^dE) = 


Fj ■ IQ-s 
cm^ s sr GeV 



(4) 


where the factors are (non-negative) adimensional co¬ 
efficients and a is the spectral index. We use the value 
a = 2, expected on theoretical basis, and find the follow¬ 
ing expressions for number of shower and track events. 

Ns = 8.4 X Fe -I- 0.9 X -I- 6.3 x Er 

Vt = 3.7 X F^ ( 5 ) 

The track-to-shower ratio is then, 


^ ^ _ 

Ns 2.3-2.0e;.-0.6er 


( 6 ) 


where we introduced the flavor fractions at Earth (i.e., 
in the detection point), defined as: 


fe = Fe/Ftot 


(7) 


with Ftot = Fe + F^ + Fr, and we considered that = 
1 — — ^ 7 -. The numerical coefficients of eq. § depend 

mildly on the spectral index, as quantified later. 

The effect of neutrino oscillations. For neutrinos 
travelling over cosmic distances, the simplest regime 
(Gribov-Pontecorvo’s [9]) applies and the oscillation 
probabilities Pm' are energy independent. The flavor 
fractions at Earth are thus given by 

= 

t i=l,3 


e? = F°/Ftot ( 8 ) 

where Ff indicates the adimensional flux normalizations 
before oscillations-see eqs. ( |4|7| ). It is generally expected, 
see e.g. [Mia that a cosmic population is characterized 
by a flavor content (^e : £,tj. ■ Ct) ~ (1/3 : 1/3 : 1/3) 
independently on the specific production mechanism. In 
this case, the track to-shower ratio in IceCube is, 

Nt 

fi = 0-24 ( 9 ) 

as can be calculated from eq. ([^ . If we consider a spectral 
index a ^ 2 , this prediction is only marginally affected 
being approximately Nt/Ns = 0.24 -|- 0.08 {2 — a). 

The equipartition of neutrino flavors at Earth is, how¬ 
ever, only an approximation which is no longer adequate 
after IceGube data: A certain imprint of the neutrino 
production mechanism does remain. It is important to 
exploit the track-to-shower ratio observed by IceCube to 
discriminate neutrino origin. To explore this possibility 
on realistic grounds, it is necessary to quantify the rel¬ 
evance of uncertainties in oscillation parameters for the 
predictions of Nt/Ns. We note that the probabilities 
Pw have a non-linear dependence on the neutrino oscil¬ 
lation parameters and, as a consequence, the errors in 6 * 12 , 
6 * 13 , 6*23 and 6 cannot be propagated linearly. Moreover, 
the allowed regions for 023 and S parameters have com¬ 
plicated structures that cannot be correctly described by 
assuming gaussian dispersions with the quoted la errors. 
We overcame these difficulties by constructing likelihood 
distributions of sin^ 0 i 2 , sin^ 0 i 3 , sin^ 023 and S from the 
N-yfl profiles given by [13]. Namely, we assume that the 
probability distributions of each parameter are provided 
by £ = exp(—Ax^/2). Then, we combine the various 
likelihood functions assuming negligible correlations and 
we determine the probability distributions of N^/Ns by 
MonteGarlo extraction of the oscillation parameters. We 
consider four specific assumptions for the flavor compo¬ 
sition at the source (^g ^ ^ C?) which are relevant for 

the interpretation of observational data because related 
to specihc production mechanisms. We consider 

i) (1/3 : 2/3 : 0) for tt decay (yellow); 

ii) ( 1/2 : 1/2 : 0 ) for charmed mesons decay (blue); 

*£/ (1 : 0 : 0) for /3 decay of neutrons (green); 

iv) (0:1:0) for tt decay with damped muons (red), 
where we made reference to the color code used in Fig. [l] 

Fig. [3 summarizes our results. The left panel is ob¬ 
tained by using the distribution of the oscillation param¬ 
eters corresponding to the assumption of normal hierar¬ 
chy (NH), while the right panel corresponds to the case of 
inverse hierarchy (IH). We see that Nt/Ns distributions 
are well separated when different neutrino production 
mechanisms are considered. This means that a precise 
determination of Nt /Ns could provide hints on the neu¬ 
trino origin, even with the present knowledge of neutrino 
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mixing parameters. From the neutrino physics point of 
view, large contributions to N^/Ng dispersions are pro¬ 
vided by the S and 023 parameters. Finally, our results 
indicates that the flavor composition of cosmic neutrinos 
cannot be used to learn about neutrinos, unless the neu¬ 
trino production mechanism is independently identified. 

For the purposes of our discussion, it is finally impor¬ 
tant to note that the track-to-shower ratio has a limited 
range of possible values, if neutrinos have cosmic origin. 
If we take the best fit oscillation parameters and assume 
a spectral index a = 2 , we obtain 


0.15 < —— < 0.30 [expected from cosmic origin] 

Ns 

( 10 ) 

The minimal value, obtained for neutron-decay (i.e. = 

= 0 and = 1 ), matters for the claims of a possible 
muon deficit problem in IceCube. If we vary the spec¬ 
tral index, this interval shifts by ~ =F 10%. The oscil¬ 
lation parameters affect slightly these expectations; e.g. 
for the lowest (resp. highest) value of sin^023 = 0.385 
(resp. =0.644) [13], the interval of Eq. (10) narrows to 
[0.16,0.27] (resp., widens to [0.09,0.43]). 


III. DATA ANALYSIS 

General considerations. In the energy window 
60TeV < Edep < 3PeV, 20 events have been observed, 
consisting of 16 showers and 4 tracks events, against an 
expected background of ~ 3 events from atmospheric 
muons and neutrinos. By performing a likelihood 
fit, an isotropic astrophysical component with E~^ 
spectrum and flavor composition (1/3 : 1/3 : 1/3), 
as expected due to neutrino flavor oscillations (see 
e.g. HDHUj)) is extracted at 5.7 ct confidence level [Tj. 
Namely, the best fit astrophysical neutrino flux is given 
by E^<^>e{E) = (0.95 ± 0.3) x 10"® GeVcm-^ s"! sr"!, 
where the index i = e, t refers to the neutrino flavor. 

New data and analyses confirm the evidence for a cos¬ 
mic neutrino population. Recently, a new technique was 
developed that permits to isolate events starting in the 
IceCube detector down to ~ ITeV and to observe as¬ 
trophysical neutrinos (in the southern sky) with ener¬ 
gies as low as 10 TeV |5]. Even more interesting, an 
independent analysis of the spectrum of muon neutri¬ 
nos passing through the Earth has confirmed the exis¬ 
tence of an astrophysical component. Analyzing the same 
period of the HESE analysis, an excess of high energy 
muon tracks is observed, that was fitted by assuming an 
astrophysical muon neutrino flux equal to E'^ <i>^(i?) = 
(1.01 ±0.35) X 10-®GeVcm-2s-isr-i [T1[TS|. 

The track-to-shower ratio. The set of events observed 
by IceCube in three years of data taking between 60 
TeV and 3 PeV consists of a total number of tt-t = 4 
tracks and ns = 16 showers. These includes on average 

= 2.1 and 63 = 0.7 background events expected from 
atmospheric neutrinos (1.7 tracks and 0.7 showers) and 


muons (0.4 tracks and no showers) [T|. In the above esti¬ 
mates, we assume that the prompt atmospheric neutrinos 
give negligible contributions, as it required by the spec¬ 
tral and arrival angles distributions of IceCube events. 
The number of tracks iVx and showers Ns which have to 
be ascribed to cosmic sources can be estimated from the 
Poisson likelihood functions: C{Ni) cx A± x e~^' where 
Ai = Vi ± bi and the index i = T, S is used to refer to 
track and shower events. By using the above data, we 
obtain Vt = 3.1 ± 2.1 and Ns = 16.3 ± 4.1. Marginaliz¬ 
ing over the total number of events, we reconstruct the 
track-to-shower ratio of cosmic neutrino obtaining 

^=0.nt° f, [HESE only] (11) 


where the error was obtained by integrating out sym¬ 
metrically (1 — CL)/2 on both sides of the Nt/Ns dis¬ 
tribution using a confidence level CL = 68.3%. The 
above result can be compared with the range given in 
eq. (10) and shows that IceCube results do not contra¬ 
dict the assumption of a cosmic neutrino population. 
The large error in the reconstructed N^/Ns is due to 
the total number of tracks which is too low to drive any 
conclusions about neutrino origin. Luckily, a completely 
equivalent and independent information can be obtained 
by the recently released IceCube data on passing muons 
M- About 12 events with visible energy above 60 TeV 
have been observed which cannot be explained by atmo¬ 
spheric neutrinos and muons. In the assumption of E~^ 
neutrino spectrum, this corresponds to a flux normal¬ 
ization F., = I.OI ± 0.35 that can be translated into a 
number of tracks from cosmic neutrinos by using eq. (|m . 
We obtain Nt = 3.7 ± 1.3 which is perfectly compatible 
with the value Nt = 3.1 ± 2.1 obtained from the HESE 
analysis, but is affected by factor ~ 2 smaller error. We 
include also this information in our analysis by construct¬ 
ing a combined likelihood, given by the product of the 2 
Poisson likelihoods for Nt and Ns and of the Gaussian 
likelihood for F/. We then extract the bound: 


Nt 

=0.18 

Ns 


+0.13 

-0.05 


[all data] 


( 12 ) 


by taking into account the equivalence between F/ and 
Nt expressed by eq. ([^ and marginalizing with respect to 
the total number of events. The likelyhood distribution 
for the track-to-shower ratio of cosmic neutrinos is shown 
by the shaded region in Figjg 


IV. DISCUSSION AND SUMMARY 

Figj^shows clearly that: i) there is no tension between 
the present observational results and the assumption of 
a cosmic neutrino population, being the central observa¬ 
tional value in the middle of the expected region; ii) there 
is no clear preference for a specific neutrino production 
mechanism, being the observational error comparable to 
the difference between the various predictions. 
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Normal Hierarchy Inverted Hierarchy 




FIG. 1. Expected track-to-shower ratio of cosmic neutrinos for the four production mechanisms described in the text. The 
distributions show the effect of uncertainties in the neutrino oscillation parameters. The left (resp. right) panel is obtained for 
normal (resp. inverse) hierarchy. The shaded region is the likelihood corresponding to Eq. {1^. 


Our results are substantially different from those ob¬ 
tained by mu]. This is partly due to the inclusion of 
the data on passing muons m, partly to the fact that 
[1[5] include in their analysis the HESE IceCube data 
between 30 TeV and 60 TeV. Following IceCube, we do 
not consider this region which is background dominated 
and much less valuable to extract the signal. 

Below 60 TeV, IceCube observes 16 events, consisting 
of 4 tracks and 12 showers [1]. The sum of tracks and 
showers agrees with the expectations but there is a deficit 
of track events (the uncertainty on the background muon 
rate is, however, 50%) and an excess of shower events. If 
one follows [HE] and supposes that most of the 12 show¬ 
ers are due to cosmic neutrinos, then Ns > 50 shower 
events are expected above 60 TeV [T], in severe contrast 
with the observational results. In other words, this posi¬ 
tion is untenable if the spectral distribution of the events, 
not discussed in [3], is considered. One possible expla¬ 
nation of the track deficit at low energy is that few 


CC interactions were erroneously identified as showers 
since the muon track was missed (e.g., for events occur¬ 
ring close to the boundary of the fiducial volume). Our 
results, expressed by Eq. (12), are stable with respect 
to a possible track misidentification systematical error. 
Indeed, above 60 TeV, the number of expected show¬ 
ers is much larger than the rate of CC interactions 
(and thus the erroneously identified events have a small 
relative importance on Ns). Moreover, Vt is well es¬ 
timated by passing muon data [14] which are free from 
track misidentification systematics. 

To summarize, the HESE events observed by IceCube 
above 60 TeV are consistent with the hypothesis that cos¬ 
mic neutrinos have been seen. The same is true for pass¬ 
ing muon events El- The flux of the cosmic muon neu¬ 
trinos can be determined reasonably well. The analysis 
of the present data gives a track-to-shower ratio, eq. (121, 
that agrees with that expected for a cosmic population, 
eq. (101. The initial neutrino flavor cannot be yet probed: 
indeed, all production mechanisms are allowed. 
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